Intershell resistance in multiwall carbon nanotubes: A Coulomb drag study 
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We calculate the intershell resistance R21 in a multiwall carbon nanotube as a function of tem- 
perature T and Fermi level £p (e.g. a gate voltage), varying the chirality of the inner and outer 
tubes. This is done in a so-called Coulomb drag setup, where a current /i in one shell induces a 
voltage drop V2 in another shell by the screened Coulomb interaction between the shells neglecting 
the intershell tunnelling. We provide benchmark results for R21 ~ V2/I1 within the Fermi liquid 
theory using Boltzmann equations. The band structure gives rise to strongly chirality dependent 
suppression effects for the Coulomb drag between different tubes due to selection rules combined 
with mismatching of wave vector and crystal angular momentum conservation near the Fermi level. 
This gives rise to orders of magnitude changes in R21 and even the sign of R21 can change depending 
on the chirality of the inner and outer tube and misalignment of inner and outer tube Fermi levels. 
However for any tube combination, we predict a dip (or peak) in R21 as a function of gate voltage, 
since R21 vanishes at the electron-hole symmetry point. As a byproduct, we classified all metal- 
lic tubes into either zigzag-like or armchair-like, which have two different non-zero crystal angular 
momenta rUa, m;, and only zero angular momentum, respectively. 



I. INTRODUCTION 

A. General considerations on nanotubes 

Carbon nanotubes are widely recognized as being 
among the most promising materials for future nanotech- 
nology applications. Furthermore, they are of fundamen- 
tal scientific interest due to several unique electronic, me- 
chanical and thermal properties. - These properties often 
depend on the microscopic details of their composition, 
e.g. the way the graphene sheets are rolled into tubes and 
whether one has a single or multiwall carbon nanotube or 
a rope or bundle of these. Electrical transport measure- 
ments have shown a tendency for ballistic transport in in- 
dividual singlewall carbon nanotubesSi2i4 (SWCNT) and 
diffusive transport in multiwall carbon nanotubes^'S'fi'i 
(MWCNT), but this issue is not completely settled yefe^ 
and seems to depend crucially on the contacts to the 
tubes and the amount of defects and impurities in and 
near the tube. Many experimentsS'ifl'ii'iSiiSiiiiii have 
explored the Coulomb blockade regime, where the tube 
can be treated as a quantum dot, due to poor elec- 
tric contact. More recently, better electrical contacts 
have been achieved which gives larger conduc- 
tance, approaching the predicted and a coherent (or 
Landauer-Biittiker-like) regime is thereby reached. Pal- 
ladium seems to be a promising candidate for good fu- 
ture ohmic contacts Another interesting feature of 
carbon nanotubes is their one-dimensional nature, which 
may have profound consequences on the basic physical 
phenomenology for their description: SWCNT's have 
been predicted to be Luttinger liquid o^^i^" and some ex- 
perimental evidence exists^liS^ even though other inter- 
pretations have been suggestediSi Whether MWCNT's 
are Fermi or Luttinger Liquids has been investigated ex- 
tensively experimentallji2iiSSiS£ and theoreticallji^i and 
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FIG. 1: (Left): The experimental setup to directly measure the 
Coulomb drag effect in a MWCNT. The intershell resistance is 
R21 = V2//1. (Right): The basic mechanism in the intershell 
resistance in a drag configuration: the intershell e — e interaction 
and thereby momentum transfer to induce the voltage drop V2. 



seems to depend on the situation, but is still subject to 
debate. Also in ropes the situation is not clear yeti^S 



The structure of this paper is as follows. We be- 
gin by introducing the intershell resistance problem in 
MWCNT's and our approach to it in section ITbI In sec- 
tion II ("I we review the basic qualitative features of our 
theory of the intershell resistance using a Coulomb drag 
setup. Sections ^ and IIIII are devoted to a summary 
of the band structure and a calculation of the screened 
Coulomb matrix element including the important sup- 
pression rules for backscattering in metallic tubes, and 
in section |^ we indicate how the standard transresis- 
tance formulae are modified in the nanotube configura- 
tion. Sections FVl I VII and I VIII give our results for several 
different nanotube combinations. Details of the nanotube 
band structure and the screening model including the 
band structure are found in the appendices. 
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B. Intershell resistance in MWCNT's 



Let us now consider electron interaction and trans- 
port in the concentric tubes (or shells) in a MWCNT. 
Yoon et al^ have argued theoretically that the in- 
tershell tunneUing of electrons is vanishingly small be- 
tween both commensurable and incommensurable long 
defect-free MWCNT's. Lack of intertube tunnelling 
is also expected in nanotube ropesi^tt Furthermore, 
Aharonov-Bohm experiments^ indicate that current only 
flows in the outer tube in a MWCNT. Another ex- 
periment by Collins et alM- supports this picture and 
finds no leakage between the shells in the low-bias 
limit. This is concluded by removing the shells in a 
MWCNT one by one and measuring the gate voltage 
response of the remaining MWCNT after each shell re- 
moval. Other shell removing experiments has also been 
reportedi^2i22i2ii2i Furthermore, Cumings et alM have 
demonstrated relative motion between the inner and 
outer shells in a MWCNT indicating that the shell are 
weakly coupled by van der Waals forces. In addition 
to Yoon et al^ also other theoretical papers have calcu- 
lated the intershell resistance using tunnelling as the only 
mechanism ii?7ii?ffn'?i?i"'fti"^i|"'^i"^i?i"'"^i"^'^ For example Roche et 
a/i^l*22i2i have considered the time evolution of a wave 
packet initially on the outer tube in a disorder-free 
MWCNT including tunnelling in a tig ht-binding approxi- 
mation. This is not in contrast to |23| due to the locaHza- 
tion of the wave packet of Refs. [37.38, 3^f22,. Using den- 
sity functional theory (DFT), Hansson et aZ.41 consider 
concentric armchair tubes and find no essential change 
in the conductance steps for a baUistic MWCNT, when 
the intershell tunnelling is turned on and off. Ref. j4^|4j|| 
also model intershell tunnelling by DFT. Very recently, 
experiments with a MWCNT with 11 contacts on the 
outer tube, where a current is driven though some of the 
tube and a voltage drop is measured elsewhere on the 
tube, have been published4S Using a transmission line 
model, information about the intershell conductance is 
deduced. 

In the present paper, we approach the intershell resis- 
tance in a MWCNT from a different point of view: We 
consider the intershell resistance R21 from the electron- 
electron (e — e) interaction between the shells neglecting 
tunnelling, i.e. in a Coulomb drag configuration. In gen- 
eral. Coulomb drag^2»^ means that moving charges in 
one subsystem (the drive subsystem) exchange momen- 
tum (and other quantum numbers) with carriers in a 
nearby subsystem (the probe or drag subsystem) thus 
exerting a drag force on the probe, inducing a current, or 
a voltage, in the probe (see Fig. QJ. Here the intershell 
or transresistance R21 = V2//1 is found as a function 
of gate voltage (i.e. Fermi level e-p) and temperature T, 
varying the chirality of the inner and outer tubes. Once 
the chiralities of the tubes are chosen, our theory has 
no remaining free parameters. Coulomb drag is a unique 



transport measurement in the sense that the R21 is domi- 
nated by the intershell Coulomb interaction4i Therefore 
serious attention to the intershell Coulomb matrix ele- 
ment and the use of proper Bloch states of the individual 
tubes is necessary. As will be seen below, the effects of 
including the band structure (and the underlying symme- 
tries of the constituent nanotubes) are absolutely crucial, 
leading to orders-of-magnitude changes in the intershell 
resistance, occasionally also reversing its sign. Further- 
more, the present work also gives a new source of friction 
against relative motion of concentric tubes, which could 
be considered in the context of using MWCNT as GHz 
nano-mechanical oscillators 

A direct measurement of the intershell resistance in a 
Coulomb drag setup (ig.GJ requires independent contacts 
on an inner and an outer tube, a difficult but possible 
technological achievement in the light of the resent shell 
removal experiments i^i^^SiS^iSiiSi As a model, we consider 
two shells, but our considerations can be extended for 
many shells. Also, a direct growth of double wall tubes 
seems feasible 1^ 

Coulomb drag has been an extremely success- 
ful tool in studying interactions in coupled quan- 
tum well a^^i^"'i^^i^^i^'''i^^ (notably in the quantum Hall 
regimesSifii) , and indeed it was realized very early that 
Coulomb drag between Luttinger liquids would be an 
important object to studySMiJ2i£2iSi£^ These stud- 
ies focused on Coulomb drag on either crossed or adja- 
cent subsystems, and used very simple models for the 
Coulomb interaction. Several interesting theoretical pre- 
dictions emerged from these papers, some of which may 
have been confirmed experimentallyiiS We work in the 
Fermi- liquid framework using Boltzmann equations. We 
think that it is important to establish a clear picture of 
what one expects within this simple model before turning 
to strongly interacting theories. Note that our approach 
also gives valuable information about drag between par- 
allel tubes. 



C. Nanotube Coulomb drag - qualitative features 
of the theory 

As explained in detail in subsequent sections, the 
transresistance or intershell resistance R21 is computed 
from the expression 

i?2ioc jY,iSRf\Vi2?A{T)F^^~>F^^\ (1) 

where the integration is taken over transferred momen- 
tum and energy in the intershell interaction, and the 
summation includes all involved bands and other quan- 
tum numbers required to specify the states. A{T) is a 
thermal factor, V12 is the screened intershell Coulomb 
interaction, and the F-functions for the two subsystems 
account for the available phase-space for electronic scat- 
tering. Of crucial importance is the factor SR account- 
ing for the selection rules (or rather suppression rules) 
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FIG. 2: The two categories of metallic nanotubes: Armchair-like 
(AL, left) and zigzag-like (ZL, right). The AL bands near = 
have zero crystal angular momentum m = and 11 = ±1, where 
fco = ^f^- The ZL tubes have doubly degenerate bands crossing 

Ep = 0, i.e. for each § = ±1 we have either tria = ^"j"™ (mod n) 
or rrii, = '^"^+^ (mod n), where ma 7^ rrv;,, (n = gcd(n, m)). The 
thin lines are the tight-binding bands near Sp = for a (n, n) tube 
(with |T| = a) and a (3m, 0) tube. 



stemming from the intershell Coulomb matrix element 
between the Bloch states. (In the final formula some SR 
is incorporated into the F-functions) . As known from 
experimentali^ and theoretica l^^i^SiTft studies, backscat- 
tering between the linear bands in metallic tubes by im- 
purities with slowly varying potentials are strongly sup- 
pressed leading to very long mean free paths. The se- 
lection rules for intershell Coulomb interaction lead to a 
similar suppression, which depends strongly on the inner 
and outer tubes' chirality. A detailed analysis of these ef- 
fects is one of the central tasks of the present article. The 
structure of Eq. Q is much richer than its counterparts' 
for coupled quantum wells due to the rather complicated 
band structure combinations of the various MWCNT's. 



II. CARBON NANOTUBE BAND STRUCTURE 

In appendixE] we give a detailed account for the band 
structure of a SWCNT with chirality (n, m), since it 
turns out to be of crucial importance to the intershell 
Coulomb matrix element and thereby also for the drag. 
Here we only outline the important points of the band 
structure used later. 

The carbon nanotube band structure can be found by 
applying periodic boundary conditions to the band struc- 
ture of a single graphite layer (graphene) . Graphene has 
two atoms in the primitive unit cell, so the tight-binding 
state (or Wannier decomposition) have two components 
with weights a and (3 (see Eq. (|A.5ll '). When applying the 
periodic boundary condition the wave vector component 
around the tube kc becomes quantized into a discrete 
values, kc — -j^^-c- However, it is important to realize 
that TLc is not the crystal angular momentum m stem- 
ming from the rotation symmetry, but only related to it 
by ric = tn (mod n). (Here C is the chiral vector and n is 
the greatest common divisor of (n,m), n = gcd(n,m)). 
This is due to the non-primitive (large) nanotube unit 
cell, when using translational symmetry instead of heli- 
cal symmetryi2ii2^ 



Linearizing the tight-binding band structure around 
the Fermi level = the states and bands for metallic 
tubes near e-p become 

— ^Hvq^t and (2) 

/ \ I ( ^ i(m-7t)-^V3(n+m) \ 

"Til 2V«^+m^+m« j, (3) 

where is the wave vector along the tube measured 
from the point, where the band crosses £p — 0,{^t — k — 
^£p=o), C = is the sign of the velocity in the band and 
<^ = ±1 describes which K<; point of graphene the linear 

band originate from. Here vq — with 70 ~ 3eV and 

a = VSoc-c {o-c-c ~ 0.142nm). The metallic states can 
thus described by (fc,^,i;). Using this, we can classify all 
metallic tubes into two categories: zigzag-like (ZL) and 
armchair-like (AL) tubes, with the following bands near 
the Fermi level (shown on Fig.[2l: 

Zigzag-like: 4m = C^^ofc, m G {ma, mf,} (4) 
Armchair-like: e]^ = —Ilhvo{\k\ ~ ka), (m = 0). (5) 

Here fc S] — j^] is the wave vector along the tube, 
ko — and T is the translational vector generating the 
translational symmetry. Note that two different tubes 
can have different |T| even though they belong to the 
same category. H = ±1 originates from ^, but does not 
give the sign of the velocity, and for a (n, n) tube 11 is the 
parity in the cylindrical anglei22i24 The linearity of the 
bands near the Fermi level is, of course, well known, but 
it is important to recognize the entirely different angu- 
lar momentum quantum numbers m that characterize the 
AL and ZL bands crossing the Fermi level. Specifically, 
for AL tubes it always holds that m = 0, while for the ZL 
tube one has tUa = ^"^"^ (mod n) or tTib = (mod n) 

(note that ma 7^ mb and nta, m;, are never zero). There 
is a one-to-one correspondence between = ±1 and the 
crystal angular momentum of the linear bands. We note 
that the most commonly studied metallic zigzag and arm- 
chair tubes, with indices (3n,0) and {n,n), are of course 
special cases of ZL and AL tubes, respectively. 

III. INTERSHELL COULOMB INTERACTION 

We next consider the Coulomb interaction between 
Bloch states \knc) for electrons in different shells in a 
MWCNT. Before calculating the Coulomb matrix ele- 
ment involving products£^ of Bloch states it is useful to 
consider the less complicated problem of the impurity 
matrix element {k' n'^V (r)\knc) . The essential assump- 
tion that we use in calculating both the impurity and 
Coulomb matrix element is that the potential is slowly 
varying on the scale of the interatomic distance flc-c- 

In 

the case of impurity scattering this is a fair assumption 
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for an impurity held on the tube by Van der Waals forces 
as is often the case i7^i77 For Coulomb interaction between 
different shells it is also a good assumption, since the 
electrons do not get close enough to experience the 1/r 
singularity. 

The impurity matrix element (k'|T^(r)|k) between the 
two component Bloch states V'k(r) Eq. IIA.5p (before ap- 
plying periodic boundary conditions) is 



R,R' 



al,akjdr^*{r ~ R')V^(r)*(r - R) 

ak'/3k ydr1'*(r-R')y(r)^'(r-R-d) 
/3ii,ak y"dr**(r - R' - d)y(r)*(r - R) 

/3iJ,/3k y"dr**(r - R' - d)V^(r)*(r - R - d) 



(6) 



By using the assumption of slow variation of V(y) we can 
take the potential outside the integrals. The first and last 
term in the square bracket become (5R',RF(R)(ak''^k + 
/3i^//?k) and the second and third term are found (includ- 
ing a sum) by summing over the nearest neighbors to 

be cx so(ak'/3kT(k') -f ak/3k'Tf"*(k)). Eq. defines 

T(k). Introducing the Fourier transform of the potential 
y(k) and the reciprocal lattice vector G we find: 



(k'|y(r)|k) = .9(k, k')-^ 5] - k + G 



(7) 



where A is the surface area and the g-factor is 



g(k, k') = ak^k' + /3k/?k' 

+ so(al^,/3kT(k') 



ak/3ii,T*(k) 



(8) 



i.e. the matrix element is essentially the plane wave result 
times a hand structure factor, which we will refer to as 
the g-factor. 

To obtain the matrix element for the screened 
Coulomb interaction y(ri,r2) (suppressing the fre- 
quency argument uj in the notation) we note that 
(k'ik'2|y(ri,r2)|kik2) = (k'2|(k'jF(ri,r2)|ki)|k2), 
where i = 1,2 labels the outer/inner tube, respectively. 
Therefore we can use the impurity potential result 
Eq. 10 to obtain: 



(k'l k'2 1 y (r 1 , r2 ) I ki k2) - .91 (ki , k'l )52 (k2 , k'2 ) 
1 



^1^2 



J2 ^(k'i-ki + Gi,k^-k2 + G2), (9) 



Gi ,G2 



G 



where we have a g-factor for each system and the screened 
potential is Fourier transformed separately in both ri and 

For a {n2,m2) tube inside a (ni, toi) tube the screened 
Coulomb matrix element is found using cylindrical coor- 
dinates r = (r, 9, z) to be 



{k[n'^^k'2n'^^\V{vi,r2)\k^n^^k^n 



1 



C2 ) 



V{k[ - ki + Gi,m[ - + nj^Ui,k2 - k2 + G2,m'2 - m2 + n2U2,ri,r2), (10) 

Gi,G2 ui,U2€Z 



where L is the length of tubes, tXi = gcd(ni,mi), Gi 
(s G Z), ri is the radius^i of tube i. 



We will also need the unscreened Coulomb matrix ele- 
ment which is a function of the interparticle distance 
|ri — r2|, i.e. a function of Z1 — Z2, ^1 — 6*2, and r2, so we 
Fourier transform in the differences zi — Z2 and Q\ — 62. 



Therefore the matrix element is: 



(fci<fc^<jy"(|ri-r2|)|fcin,^/c2n,^) 



1 

-5i('^i"ci 

/ ^ / ^ y°(fci - /ci -f Gi,mi - nti -I- n^Ui, ri, r2) 



> )52 (fc2^C2 ' ^2<2 ) 

Gl,G2 Ml,tl2eZ 

X '^fcj+fcj +fe2+'-'l+'-'2 ^™'l+'^2+'''l "1 ''^l+'^2+"2"2 ' i^^) 



Note that the ± in the states Eq. IIA.llll (the £, index for 
metalhc states Eq. (I2I3II ') is suppressed in the notation 
and that this index only appears in the g-factors in both 
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Eq. and Eq. Illlll . Here we have used the crystal 
angular momentum difference in the Fourier transforms 
instead of the ric difference, since this is the physical 
(crystal) angular momentum being transferredi2^ Note 
that we have included Umklapp scattering and that the 
unscreened interaction Eq. lfTT|l has crystal (angular) mo- 
mentum conservation. Similar matrix element were con- 
sidered by S. Uryu4^ 



A. The g-factor and backscattering in metallic 
tubes 

We now consider the (7-factors and show that they con- 
tain essential information about the electronic scattering. 
The g-factor for any {n, m) metallic tube for the scat- 
tering process (fc,^,<^) (fc',^',<;') between the metal- 
lic states Eq. ^ is found by inserting Eq. ljA.lfl|l (with 

= ^T-m) and Eq. © into Eq. ©: 



Zigzag- like : Zigzag-like : Armchair-like : 




FIG. 3: The possible backscattering processes in any metallic tube 
with a slightly raised Fermi level e-p. (Left): Backscattering in a 
zigzag-like tube without crystal angular momentum change Am = 
(i.e. <; = ?') and a small wave vector \k' — A;| ~ change, which 
is suppressed by g < 10~^ from Eq. 1141 . (Center): Backscattering 
in a zigzag-like tube with crystal angular momentum change, which 
have g ~ 1 from Eq. 1151 . Here m denotes the opposite of m in the 
set {ma,mi,}. (Right): Two types of backscattering in armchair- 
like tubes: (i) A large wave vector transfer (for ? = — ?') in between 
states with the same crystal angular momentum (m = 0) and g ~ 1 
(Eq. 1151 *) and (ii) A small wave vector transfer g ~ suppressed 
by g < 10~^. Note that the distance between the points ii^fij are 
not to scale (i.e. |^ <C ^M^) and that the armchair-like bands are 

^ hVQ d I 1 I 

connected as in Fig.|2] 



where we introduced 



-t- -I- 4mn 



, . V3.?(m^-7^2) 

2{'n? + + mn) ' ^ ' 



The g-factor in Eq. llT2|l has two terms: The first paren- 
theses is the important wave vector independent scalar- 
product of (j^) from Eq. Q and the second term is a 
wave vector dependent correction term (of first order in 
so -0.1). 

As we shall show in section IIVL only backscattering 
contributes to the Coulomb drag in metallic tubes and 
we therefore need to consider all possible backscattering 
processes (^ = — ^') in any metallic tube. 

Due to the double degeneracy of the zigzag-like bands 
Eq. Q at the Fermi level, we must consider backscatter- 
ing both with and without crystal momentum exchange 
(Fig. center and left panels, respectively). 

If 1; = <j' then Am = m' — m = and from Eq. l(T2|l we 
have 



|g(fc,(;,^;fc',^,-^)| 



^/ia\k' -k\ 



(14) 



which is of order 10 ^ or less for scattering around the 
Fermi level, i.e. for \k' — fc| ~ the g-factor is \g\ — 



So- 



< 10^2 for < o.3eV. If <r = -<:' then |Am| = 



level the g-factor squared isiSi 



ill 



n + m + 4mn 
2{ii? + m? -t- rrm) 



16 V '^^ + fn^ + Trm 



(15) 



Itria — tribl 7^ and for backscattering around the Fermi 



which is I for (n, 0), 1 for (n, n) and in between for 
all other tubes. So in a zigzag-like tube we have two 
kinds of backscattering with small crystal wave vector ex- 
change q ^ (and thereby large V{q, Am)): Either 
Am = and \g\ < 10^^ or Am 7^ and \g\ 1. 
Note that the larger the Am the smaller V{q, Am). Even 
though V{q, Am) is large the small (7-factor suppresses 
the Am — backscattering. 

Consider now armchair-like tubes where the bands 
crossing £p = all have m = 0, so the small crystal 
wave vector transfer around ±^f^ have = and there- 
fore the g factor is the same as in Eq. Ill4p . i.e. \g\ < 10^^ 
suppresses this kind of backscattering (Fig. right)). If 
we on the other hand have a large crystal wave vector 
transfer backscattering (Fig. right)), then <; = — and 
the g-factor of order 1 from Eq. Ijl5|l is used. So the 
large crystal wave vector backscattering is most impor- 
tant, since the Fourier transform does not grow enough 
to compensate for the small g-factor. 

Ando et alm^t^ have used the k • p approximation 
to consider backscattering (from impurities) in metaUic 
tubes and found a result similar to Eq. but without 
the So term. The small wave vector transfer backscat- 
tering was found to be small in these papers. Klessa^ii 
has found similar results for scattering in nanotubes, 
see also for some experimental evidence of lack of 
backscattering in metallic tubes compared to semicon- 
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ducting ones. 



ing the carbon nanotube band structure with the result: 



B. Screening effects using RPA 

In appendix B, we derive the screened Coulomb inter- 
action in the random phase approximation (RPA) includ- 



J 



2ttL 



Gi ,G2 "1 ,"2 



E 



V°{k[ - ki + Gi.rrii - +niUi,ri,r2) 
ei2{k[ - ki,m[ -~m^,uj) 



r 



where ei2((?, Am, w) is the dielectric function disregard- 
ing the Umklapp processes (see Eq. IIB.lflll V Note that 
the effective non- interacting polarization Xes i ^) 
Eq. IIB.8|l entering the dielectric function contains the g- 
factors. For armchair-like tubes x^^j ^(g, Am, w) is given 



explicitly in Eq. llBJ2ll and Eq. llBJ3ll . The 

Coulomb interaction for cylindrical geometry is 

^2 



bare 



l/°(g,Am,r, 



— lAm('7''i)KA,i 
£0 



(qrj) n <rj, (17) 



where lAm{x) (KAm(a^)) is the modified Bessel's func- 
tions of the first (second) kind of order Am and eo is the 
vacuum permittivity. 



IV. THE TRANSRESISTANCE MODEL 

The transresistance R21 is now found for diffu- 
sive nanotubes using two coupled BoUzmann equations 
(i.e. Fermi liquid theory) in linear response to the applied 
electric field Ei and for weak coupling between the tubes. 
The derivation is a generalization of ref. 
(used to study bilayer systems) to the case of several 
general bands. We only sketch the derivation and the 
details can be found in Chap.3 of Ref. In order to 

simplify the notation we use ly as a collection of band 
indices for the tube. A similar formula of R21 can also 
be found using the Kubo formula and doing perturbation 
theory to second order in the intertube interaction (the 
first order DC contribution is zero)i£^ 

The coupled linearized Boltzmann equations for the 
non-equilibrium distribution functions fi(ki,Vi) (j = 1,2 
see Fig. m are: 



eiE^df{eu,u,) _ /i(fci,^i)-/"(£fc,.J 
h dk^ 

e2E2 9/"(£fe,.J 



(18) 



/2(fc2,^^2)-/°(£fe,.J 



+ 5[A,/2 = /°](fc2,i^2) (19) 



where a simple relaxation time approximation is used for 
the impurity scatteringS^, is the carrier charge in sub- 
system i and S[fi, A — /°] is the linearized collision inte- 
gral coupling the two subsystems/tubes. The assumption 
of weak intertube interaction and small external electric 
field El were used to linearize the equations and to only 
include the lowest order terms and therefore not have 
a collision integral on Eq. III8II . The linearized collision 
integral is (using the H -theorem^): 

S[h,f2 = f]{k2,^2)= (20) 

-EE E -(l'2';12)/°(e.^.J/°(£.,.J 

X (1 - /°(£A:;.;)) (1 - f{ek'2u'2)) [Mku^i) " MKA)] 

where the deviation from equilibrium ipi{k, v) was defined 
though /,(fc, v) - f(ek.) = f°{sk,){l - fieku))Mk, y) 
and w(l'2'; 12) is the transition rate for electron-electron 
scattering between the tubes found from the Fer- 
mis golden rule w(l'2'; 12) = xK^i''i^2i^2l^i2(|ri - 
v-i\)\k^v^k^v^\^8{ek^y^ + Sfc^..^ - ey^v^ - ^fe^,.^) using the 
matrix element in Eq. ifTfiji . To derive the transresistance 
i?2i = we use the coupled Boltzmann equations Ijl8|l 
and Ijl9|l with 112011 and that I2 — 0, since a voltmeter is 
placed on subsystem 2i2i After some algebraSi we get: 



-R21 — 



7reie2nin2fcgT (27r)^rir2 
1 



(27r) 



=dg 
27r 



V12 (g, , K > ^1*2 (<? + Gi , , , w) 



sinh^ ( 



2k„T) 



X F^'^\(q,u;)F^^\(q,uj), 



(21) 



where rii is the carrier density, Vi2{q,h'i,v'i,uj) = 

KXaZT '^^^^2^2) are the se- 

lection rules for the band indices such as crystal angular 
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momentum and / or parity (for armchair tubes) conserva- 
tion and i^j^'|,,(q, w) is the available (g,tj)-phase space for 
scattering in the z:te tube given by 



5D 



i ' i 



(22) 



sign(ufe^i., - Vk^+q^0 

X [f{ek^.J-f{ek^+,.0^ \g,{k,i.„k, + q,,',)\', 

where the kg are the solutions to Eki^. — £k+qv' — fiu! = 
in the FBZ of subsystem i, Vkv — ji^^^ is the velocity, 
sign(a;) gives the sign of a; (if a; = then sign(a;) = 0) and 
fi^j. is the transport mobility, which is a single subsystem 
property. Note that the F-function is periodic and odd 
in q. 

Having stated this formula a few comments and inter- 
pretations are in order. 

Firstly, we note that only backs cattering processes con- 
tribute to the drag between metallic tubes in the linearized 
band models Eq. Q and Eq. 10 , since we only have two 

velocities ±uo — ± in the metallic bands and there- 

fore the sign-function of the velocity difference before 
and after the scattering event in the i^-function Eq. Ij22|l 
makes only backscattering (i.e. Vk^„i = —Vk^+qu') con- 
tribute to the F-function. In section IIII Al we there- 
fore analyzed the g-factors for all possible backscatter- 
ing processes in metallic tubes. The interaction and 
sinh~^ ( 2k^T ) ^""^ decreasing functions of q and uj, re- 
spectively, so the importance of the phase space (i.e. the 
F-functions) in the integral decreases from the origin. It 
is worth to note that the forward scattering contribu- 
tion which for quadratic dispersion relation dominates at 
higher temperaturesjr^ thus plays no role here. If we in- 
cluded a curvature of the dispersion relation for the nan- 
otubes, we would get a correction to the results presented 
here. However, there is one subtlety hidden in this, be- 
cause if we consider Coulomb drag between short tubes, 
where the distribution functions are not relaxed to the 
Galilean invariant form assumed in but is instead 
given by a two-step distribution function, the forward 
scattering does not contribute to the Coulomb drag as 
shown in 

Secondly, we have used a quantum number indepen- 
dent impurity relaxation time Ti for each subsystem in 
Eq. I^Hl and Eq. ltT!?)l . The mobility can be shown 

to be proportional to n, i.e. fi^^ cx Ti, from a single sub- 
system Boltzmann equation (like Eq. lllSll l. Therefore 
the F-function Eq. (I22|l is Ti independent, so the trans- 
resistance i?2i is independent of the impurity relaxation 
times. So in the quasi-ballistic regime for large n the 
transresistance is still formally correct. However, there 
has been some work on drag between ballistic ID systems 
with free electron like bands using Boltzmann equations, 
where almost identical transresistance formula is foundi^ 
As a last comment, we note that Umklapp scattering 
is only possible if the tubes are commensurable due to 
the Sgi,G2 function in Eq. II21|I as also found in ref. [s^ . 
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FIG. 4: The transresistance per length (in U/fim) as a func- 
tion of the Fermi level £p (in eV) (e.g. a gate voltage). The tem- 
perature is: T = 80K (dotted), T = 150K (dashed) and T = 300K 
(full line). The dip in R21 at £p = reflects the electron- hole 
symmetry at this point. (Inset): A sketch of the situation for mis- 
aligned Fermi levels (see text). 



V. ELECTRON-HOLE SYMMETRY AND 
COULOMB DRAG 

All nanotubes have an inherited electron-hole symme- 
try from the graphene band structure for gp = 0, which 
intuitively means that there are as many electrons as 
holes for Ep = (for the precise definition see for a 
recent measurement of electron-hole symmetry see ^1,3]) • 
So there will be an equal amount of momentum trans- 
fer to (from) the electrons and holes and therefore no 
voltage difference will arise, i.e. R21 = 0, if one of the 
subsystems has electron-hole symmetry. Formally, the 
F-function can be seen to vanish at electron-hole sym- 
metry by using /°^(e) = 1 — /°(— £) (after doing the 
sum over the band indices), where fj. is the chemical po- 
tential. This has also been used to show how R21 can 
change signi^ 



Therefore, we predict that by varying the Fermi lev- 
els (either by gate voltage or doping) a dip (or peak) in 
R21 will appear due to the electron-hole symmetry for all 
kinds of tube combinations. For two concentric armchair 
tubes ((5, 5) in (10, 10)) the transresistance as a function 
of £p is shown on Fig. 31 (the range of Ep is chosen to cor- 
respond to typical experiment8Si££*£2*2SiSii2^) . Note that 
we use the simplification of having the same Fermi level 
in the two tubes in the numerical calculation, but the sit- 
uation for misaligned Fermi levels is sketched in the inset 
of Fig. 31 As indicated in the inset, we have the following 
scenario for increasing gate voltage for Ep' ^ Ep^: First 
hole-hole (h — h) scattering, then one subsystem passes 
electron-hole symmetry, i.e. R21 — 0, afterwards e — h 
scattering until the other subsystem also passes though 
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the electron-hole symmetry point. The details of the cal- 
culation are given below in section FVI Al 



VI. COULOMB DRAG BETWEEN METALLIC 
TUBES 

A. Drag between (real) armchair tubes 

Let us begin by calculating the transresistance Eq. lf2T|l 
between two concentric real (i.e. {n,n)) armchair nan- 
otubes, which have |T| — a independent of n. The band 
index f is in this case the index 11 = ±1 from Eq. J^. 
To find the F^^,{q,u!) functions Eq. lf22|l we need the 
Efc+g — huj = with the bands Eq. Q 



solutions of £^ 
and remembering that £^ should be made periodic 
by hand (in order to find two solutions and not only one). 
The sign function only gives backscattering, which is ex- 
pressed by step functions. For intraband backscatter- 
ing n' = n we have g ~ 1 (Eq. Ijl5|l ) and for inter- 
band backscattering 11' — —11 we have \g\^ — So^^^ff^ 
(Eq. itTljl ) as found in section IIII Al Therefore the F- 
functions are^S for < g < 



F['l{q,uj) 



ITIT- 

^d;'>9{-UJ + Voq) 



(23) 



- £2 



£i 



-kohvo 



-kohvo 



with El — ^{u} + voq-2vako) and £2 = |(^-wo9 + ^'ofco), 



F^l{q,uj) = -C^t 



-uj + voq) 



(24) 



£2 



£1 



Hvoko 
hvoko 



with ii — |(w + i;og — fo^o) and £2 = K"^ — wo9 + 2wofco), 



F«(<Z,^) = -C 



16 



■Q(uj - Voq + voko) 



(25) 



X 6{~lli — Voq + 2voko) 
{f{e',)-f{-e',)) + {f{e',)-f{-e',)) 



and 



[i) 2%^ 



16 



{-LO - voq + voko) (26) 



(/°(£'i)-/"(-£'2)) + (r(£:i)-r(-£'i)) 



where e'l — ^{uj + voq) and £2 = f (w — voq) and we have 
calculated the common single subsystem prefactor 



C 



2£p + fiwojT-J 
2ihvo)^ ■ 



(27) 





FIG. 5: Contour plot of the F functions for the intraband scatter- 
ing for < g < £p > and the temperature T = O.lTp. Note 
the smearing by the Fermi functions due to the temperature on 
some edges and the sharp edge at oj = voq from the step function 
e{-uj + voq). 



It is important to note that the interband -F-functions, 
F^ and F j. , are heavily suppressed compared to the in- 
traband i^-functions (shown in Fig. 01 by jgp = 5§ '^'•"g-' 
of order < 10~^ for backscattering around the Fermi 
level. Therefore, including the tight-binding states in 
the Coulomb matrix element and not just in the avail- 
able phase space for scattering as in Ref. is a very 
important effect. 

In real armchair tubes the 11 index is a parity index in 
the cylindrical coordinata£2iii and therefore the Coulomb 
matrix element has the property: 

{k[ii[k'2ii'^\v{vi,v2\)\k,n,k2ii^) = 

u[u'^n,n^{k[n[k'^n'^\v{r^,r2\)\k,n,k^n^), (28) 

i.e. the product of the parity is conserved in the interac- 
tion. Since both 11 = ±1 have m = there is no angu- 
lar momentum selection rule, so the only selection rule 
J in Eq. (EU is J{n„U[,U2,U'^) = Sn^u,,n{u',, which 
reduces the number of terms by a factor of two. Since 
V{q, Am) is parity independent in Eq. 1I2TI1 . then the sum 
over band indices for |Ti| = IT2I is: 



{F++ + + + F+_)' ^ (Fi„t,.a)2 + (Fi 



(29) 



inter ) : 



which defines the inter and intraband F functions. 
(Fintcr)^ is of fourth order in soq and therefore strongly 
suppressed compared to Flntra even though Fjntor has 
phase space for smaller q and lu. Fintra(<Z,w) is shown 
on Fig. H 

We now have all the ingredients of the transresistance 
-R21: 



1 



-d, rjvM,M\^ 



R2I ^ fi' 

L Txe'^nin2k^T 27rrir2 Jo 27r Jq " sinh^ V2FT; 

X [(Fintra(9,w))' + (Fintcr(<Z,w))']. (30) 

A numerical integration yields R21 as a function of £p 
and the temperature T, shown on Figs. and respec- 
tively. The transresistance per length is of the order 
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FIG. 6: Contour plot of the 

■Pintra function. Fintra giveS 

the phase space for intraband 
scattering in (real) armchair 
tubes. -Fintra is seen for < 
g < ^ and is odd in q and 
should be repeated periodi- 
cally with 22L as a function of 
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FIG. 8: The transresistance 
per length versus radius 

(r oc n) for armchair tubes. 
The different outer and inner 
armchair tubes are: A (5, 5) 
in a (n,n) (dots), a (6,6) in a 
(ra, n) (triangles) and a (9, 9) in 
a (n,n) (stars). The radius of 
the outer tube is: r = ^^^n 
for a (n, n) tube. Here T = 



35 300K and e-n 



0.3eV is used. 



Note the logarithmic scale. 
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FIG. 7: The transresistance per length ^^j^ (in units of f!//im) 
versus temperature F (or T/Ty (left)). The curves are obtained 
from a numerical integration of en. 1,301 for a (5,5) in a (10,10) 
tube. Curves for four different Fermi levels £p (i.e. gate voltages 
or dopings) are seen: £p = 0.006eV (Fp = 69K) (left, dashed 



0.015eV (Fp = 174K) (left, full line), 



0.15eV 



line), 

(Tp = 1740K) (right, dashed line) and = 0.3eV (right, full line). 
Note the difference in magnitude between the transresistances R2i- 



a few fl/^m. R21 is seen to be linear in T for T < 0.4Tf 
as also found for free electron like bandsifi^ For higher 
temperatures the transresistance increases or decreases 
depending on the Fermi level. Numerically, we find a 
factor of 10^ difference between the contribution to R21 
from i^inter and fintra, SO wc Can conclude that the drag is 
due to the intraband hackscattering processes. The largest 
contribution to the integral is around q — ko± |^ (see 
Fig. EJ, which corresponds to Umklapp scattering pro- 
cesses around the Fermi level, e.g. fc = fco — ^ and 

A:' = -fcn + ^ so a = fc' - A: + ^ = fco + 



Note that screening induced by the substrate could 
change the magnitude of the transresistance a small 
amount, which could be modelledSfi by introducing a new 
dielectric constant k = er-eo instead of eq in Eq. ljB.15|l 
with about 1 to 2>^ For the present case, the magni- 
tude of R21 is changed < 10%, when is increased from 
1 to 3. 

The transresistance depends on the radii of the tubes 
only via the bare Coulomb interaction Eq. IIB.1511 . Fig.|H| 
shows that R21 decreases exponentially (for n < 25) when 
keeping the inner armchair tube at a fixed radius and 
increasing the outer tube radius. For parallel 2DEG's R21 
was found to depend on the separation d as5£ R21 oc d^. 



B. Drag between armchair-like tubes 

For two general armchair-like tubes, we do not have a 
parity selection rule and in general |Ti| ^ IT2I as seen 
in Tabled Therefore we have no selection rules, but all 
other terms than Fjj^pjj)^, are of higher order in (so?)^ 
and therefore small, i.e. 



n,n,n', n: 



p{2) p{2) 

^ n,n', ^ n,n' 



Fi'lFi'l + F^'lF^'l 



^(1)^(2) 



(1) p(2) 



F^'1F_ 



(31) 

as for the (real) armchair tube case Eq. Il29|l . The F^l_ 

and F^l are the same as found in section lVI Al and shown 
in Fig.0except that a is replaced by \Ti\ (but not in the 
g-factor) . 

Since |Ti| and IT2I are different (in general), it is 
harder to conserve (crystal) momentum near the Fermi 
level for the dominant hackscattering process with mo- 



mentum transfer q 
ever, for some values of 



'^0 ^ hvo 



With k^;^ = 3^. 



How- 



Ti| and IT2I it is possible to 
conserve momentum near the Fermi level, which gives 

IT I 

rise to peaks in i?2i e.g. at |^ = 
The peaks on both sides of — 



IT2I 



1 as seen in Fig. El 
1 are 



|Ti| 



IT1IIT2 



6et 



irhvo 



(32) 



corresponding to fcg^^ ± ^ 



_ 1.(2) 



T 1^ (see inset (a) 
in Fig. EJ. These peaks have i?2i < 0, since they corre- 
spond to a resonance between a electron-like and a hole- 
like hackscattering in the sense that a hole-Hke (electron- 
like) hackscattering takes place in a hole-like (electron- 
like) band with sign(wfc) = — sign(fc) (sign(wfe) = sign(fc)) 
in the FBZ. The peaks around = ^ and 2 are found 
in the same way by taking the hackscattering processes 
q ~ 2fcQ*' ± 1^ into account. If the radii of the tubes 
are different, then the magnitude of R21 will change (see 
Fig. but the signs and positions of the peaks are the 
same. The peaks are broadened by increasing tempera- 
ture and the positions of the peaks depend on £p as seen 
e.g. from Eq. ^ (except for |^ = i, 1 and 2). The 
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FIG. 9: The transresistance per length as a function of the 
ratio of the translational vectors length |Ti|/|T2| for two armchair- 
like tubes. The peaks corresponding to different scattering pro- 
cesses are seen as explained in the text. Numerically, we use 
IT2I = a, radii as for a (5,5) in an (10, 10) tube, T = 300K and 
Ep = 0.3eV. If the tubes have a different radius, only the magnitude 
of the peak is changed, see Fig.|s] (Inset (a)): The scattering pro- 
cesses in tube 1 and 2 leading to the peak at |Ti|/|T2| ~ 1.28. Note 
that the backscattering processes are electron-like and hole-like, re- 
spectively, so i?2i < 0. (Inset (b)): Peaks around |Ti|/|T2| = 1/2. 
Note the difference in scale. 



1. If Am ^ is not possible, then only Am = 
backscattering for q ~ |2z. jg present, but this is strongly 
suppressed by the small g-factor and so is the drag. So in 
this case the small wave vector transfer forward scatter- 
ing (for non-linearized bands) could become important, 
but in any case the effect is small. An example is the 
drag between two (real) metallic zigzag tubes (see table 

IB. 

2. If Am 7^ is possible, then this process is the domi- 
nant, even though there is a small suppression (compared 
to the g-factor) from having Am 7^ in the Fourier trans- 
form 1/12(9, Am, w), which is smaller the larger Am. An 
example is a (12, 15) in a (15, 18), which has an angular 
momentum exchange of Am = ±1. 

Furthermore, there are no peaks in R21 as a function of 

IT I 

{tpH as for the armchair- like tubes, since the transferred 
IT2I ' 

crystal wave vector q ~ is independent of |Ti|. 

Prom the same principles as used above, we find the 

drag between zigzag-like and armchair-like tubes to be 

strongly suppressed. 



VII. COMMENTS ON THE DRAG BETWEEN 
SEMICONDUCTING TUBES 



situation of varying |Ti| and IT2I is similar to varying 
the densities in the parallel 2D systems jSi Note that if 
we have a tube configuration corresponding to a nega- 
tive dip in Fig. |2l {R21 < 0) , then this tube configuration 
will have a peak instead of a dip as a function of the gate 
voltage. 

Summarizing, Coulomb drag between armchair-like 
tubes is strongly dependent on the magnitude of the 
translational vectors |Ti| and IT2I and can lead to both 
negative and positive transresistance. 



C. Drag between zigzag-like tubes 

Consider the drag between two zigzag-like tubes, where 
the V index is ^ = ±1 and m € {ma,mb} from Eq. I^J. 
The backscattering i^-function, F^J^ _^rn'' ^ form sim- 
ilar to -FnLn for armchair-like tubes (Eq. lf25|l and lf2fill ). 
where the important part is the backscattering around 
the small q ~ This backscattering can be both with 
(Am 0) and without (Am = 0) exchange of crystal 
angular momentum with the g-factors 

|5(Am = 0)p cx (soag)^ and |g(Am7^0)p~l (33) 

found in section UlI Al Since there is crystal angular mo- 
mentum conservationSSi it depends on the combination of 
the zigzag-like tubes (and their m^ and mf,) whether the 
Am 7^ backscattering is possible or not, so we have two 
very different cases: 



If the Fermi level for a semiconducting tube is shifted 
into the conduction (or valence) band, then the drag pro- 
cesses are within a single band (i.e. Am = 0) similar to 
a quadratic band for small tubes, where there are few 
bands with large separation. Here both the small q for- 
ward scattering and the large q backscattering processes 
will contribute to the drag. We can calculate the g factors 
in the same way as for the metallic tubes and for intra- 
band scattering they are of order one. However, the mag- 
nitude of the backscattering momentum transfer around 
the Fermi level has to be approximately the same in the 
two tubes in order to satisfy momentum conservation. In 
general, this is not the case. 

If we deal with larger tubes more bands can come 
into play and thereby more scattering possibilities ap- 
pear than captured in the single band quadratic model 
(see ref. for a discussion on scattering in larger 
MWCNT's). This is also the case of larger metallic 
tubes. Coulomb drag in the quadratic model with more 
bands (with different angular momentum along the tube) 
for tubes of semiconducting material are considered in 
ref. 113 ■ 



VIII. SUMMARY 

We have considered the intershell resistance i?2i origi- 
nating from the intershell Coulomb interaction neglecting 
tunnelling, i.e. in a Coulomb drag configuration. 

For any tube combination we predict a dip or peak 
in i?2i as a function of gate voltage, which should be 
experimentally observable. The dip (or peak) is due to 
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the electron-hole symmetry of the carbon nanotube band 
structure. Whether R21 has a dip or peak depends on the 
sign of i?2i, when both systems have Fermi levels above 
the electron-hole symmetry point. 

The order of magnitude and sign of R21 were found 
to depend crucially on the chirality and Fermi level mis- 
matching of the two tubes. The magnitude of R21 can 
reach ~ 50ri//xm under favorable circumstances. The 
origin of the drastic change in magnitude between differ- 
ent chiralities is the suppressed backscattering due to the 
Coulomb matrix element between Bloch states combined 
with the mismatching of wave vector and crystal angu- 
lar momentum conservation near the Fermi level. R21 
was found to be linear in temperature for low temper- 
atures (compared to Tp), just as for a single quadratic 
band. To facilitate the analysis, we classified all metaUic 
tubes in two categories: zigzag-like or armchair-like, and 
described their crystal angular momentum properties. 

Throughout the paper, we use Fermi liquid theory to 
describe the Coulomb drag in the MWCNT's, which gives 
a benchmark result for comparison to future experiments 
and Luttinger liquid theories of drag in MWCNT's. The 
effects considered in this paper should be helpful in inter- 
preting future measurements of the intershell resistance. 
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APPENDIX A: ENERGY BAND STRUCTURE 
OF THE CARBON NANOTUBES 

We will now give a rather detailed discussion of the 
band structure of carbon nanotubes, since the intershell 
Coulomb interaction matrix element turns out to depend 
critically on the Bloch states of the two tubes due to the 
two atom primitive unit cell (of a graphite layer) as seen 
in section Unl 

The carbon nanotube lattice can be thought of as a 
wrapping (i.e. a conformal mapping) of a graphite layer 
into a tube. The wrapping is preformed such that the 
chiral vector C = nai -I- ma2 becomes the circumferen- 
tial of the (n, m) nanotube and this determines the lattice 
completely32a (Here ai = f (\/3, -1) and a2 = f (\/3, 1) 
are graphene lattice vectors and a = |aj| = VSac-c, 
where Gc-c is the inter atomic distance). 



Any (n, m) nanotube has three symmetries: A dis- 
crete translational symmetry along the tube, a discrete 
rotational symmetry around the tube axis and a heli- 
cal symmetry (i.e. a screw operation). These symmetries 
give rise to the three corresponding quantum numbers: k 
(crystal wave vector along the tube), m (the crystal an- 
gular momentum component along the tube) and k (he- 
lical quantum number). Only two of these symmetries 
(quantum numbers) are needed to label the eigenstates, 
since the symmetries are not independent Convention- 
ally translational symmetry is used to label the states, 
but this does not use the smallest possible unit cell and 
can therefore give many bands in the first Brillouin zone 
(FBZ) with the same angular momentum. 

Any carbon nanotube can be generated from a primi- 
tive two atom unit cell using only discrete rotations and 
discrete screw operations and thereby giving (general- 
ized) Bloch states |Ktn)iii*2^ The advantage of using this 
method is that each energy band (as a function of k) has 
its own crystal angular momentum m. The discrete rota- 
tional symmetry is generated by the vector Cn along C 
giving the smallest possible rotation leaving the lattice 
invariant, i.e. 

Ti vn 

Cn = — ai H a2, where n = gcd(n,TO), (A.l) 

i.e. n is the greatest common divisor of n and m. So 
a given (n, m) tube has crystal angular momentum m £ 
{0, 1, ... , n— 1}. The disadvantage of using the symmetry 
adapted Bloch states |Km) is that n is in the direction the 
generator H of the helical symmetry, which is general is 
different for different chiral vectors. 

If we instead use the (often much) larger translational 
unit cell the states can be labelled by k s] — ^,fr], 
where T generates the translational symmetry (the trans- 
lational vector) and is given bjiiffi 



(2m -|- n)ai — (2n + m)a2 
gcd(2m -\-n,2n + m) 



(A.2) 



Since we do not use the primitive unit cell in this case, 
but a larger translational unit cell, we get a smaller FBZ 
and thereby more bands in the FBZ than there are crystal 
angular momentum quantum numbers. 

The conventional way to obtain the band structure for 
a isolated singlewall (n, m) nanotube using the transla- 
tional unit cell is to apply periodic boundary conditions 
on the two dimensional graphene tight-binding stateiSi 
'0k(r) along the circumferential C of the tube22ifi£, i.e. 



V'k(r + C) 



ikC 



V'k(r) = i^kir) ^ k - C = 27rnc, (A.3) 



where ric is an integer in {0, 1,2, ... ,J\f — 1} with TV = 
2(n +m +nm) ^ ^ being the number of (two atomic) 



gcd(2m+n,2n+m) 

graphene unit cells in a translational unit celliiffi Thereby 
the Uc labels the bands (as a function of k) using the 
translational unit cell. One disadvantage of using this 
larger translational unit cell is, that is not the crys- 
tal angular momentum, but only related to the actual 
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Primitive unit cell Zone Folding procedure Translational unit cell 




FIG. 10: (Left): The m = 1 band for a (5, 5) tube in the FBZ of 
the primitive unit cell as a function of k G] — it,tt]. (Center): The 
m = 1 band is pushed into the smaller FBZ of the translational 
unit cell by using k = ^ + and that nc = 1 and ric = 6, since 
m= 1. Note that the band is symmetrical around ^, since m= 1. 
(Right): The band structure for the translational unit cell. Both 
bands have crystal angular momentum m = 1, but indices ric = 1 
and ric = 6. 



physical crystal angular momentum m by: 

ric — m (mod n). (A. 4) 

Furthermore, we can connect the description of the band 
structure using the primitive unit cell and the transla- 
tional unit cell by k = k • H, i.e. k depends on both k 
and ric^ An example is given in Fig. 1101 

To do a tight-binding calculation for grapheneiiii it is 
essential that the unit cell of graphene has two atoms, 
so the tight-binding state (Wannier decomposition of the 
eigenstate) has two components^S^: 

Mr) ^^T. ^^■''(ak^lr - R) + /3k*(r - R - d)), 

(A.5) 

where VE* is a 2pz orbital (Wannier function) localized at 
each atom, R = niai -I- n2a2 (ni,n2 € Z) are lattice 
vectors, N is the number of unit cells in the layer, d = 
i (ai -|-a2) is the vector between the two atoms in the unit 
cell and a^, /3k are functions to be determined by the 
tight-binding calculation. To find the energy we insert 
ipk{v) in H%l)\^(v) = ekV'k(r) and obtain a 2 x 2 matrix 
equation: 

(S;S:)(s)-'(S;£)(s)' <"-^' 

where Hij, Sij are the overlap integrals with and with- 
out the Hamiltonian found in the nearest neighbor tight- 
binding approximation to be: 

Hi2 - -7o(l + e-*-"^ + e-*-^) = -7oT(k), (A.7) 
<S'i2 = soT(k), Su = 1 and Hu = eo, (A. 8) 

where the value of the overlap integral is 70 — 3eV, the 
overlap of the orbitals areiffi sq ~ 0.1 and £0 is the energy 
of the orbital, which is set to zero. Here the convention 
70, So > is used. In the literature a range of different 
values is used (e.g. 70 ~ 2.5 - 3.i),i"3,in4,in.5,in6 
diagonalization of Eq. IIA.fijl we find (for T(k) ^ 0): 



Chirality 


Type 


ma 




n 


\T\/a 


D/a 


(n,n) 


AL 








n 


1 


TT 


(7,4) 


AL 








1 




793 

TT 


(15,6) 


AL 








3 




3V^ 

TT 


(8,23) 


AL 








1 


^/259 


TT 

5V39 

TT 


(10,25) 


AL 








5 


\/T3 


(n,0) 


ZL 


2n 
3 


n 

■3 


n 


V3 


n 

TT 


(9,6) 


ZL 


2 


1 


3 


V57 


3^19 

TT 
TT 


(6,21) 


ZL 


2 


1 


3 


V20T 


(18,12) 


ZL 


4 


2 


6 


VEr 


TT 
I2V7 


(12,24) 


ZL 


4 


8 


12 


V21 



TABLE I: Examples of armchair-like (AL) and zigzag-like 
(ZL) metallic tubes, i.e. all kinds of metallic tubes. For the AL 
tubes the difference in the length of the translational vector T 
and the diameter D is seen and for the ZL tubes we note the 
variety of the crystal angular momentum rria — ^"3 (mod n) 
and rrii, = ^"y"" (mod n) of the bands crossing the Fermi level 
(ep — 0). Numerically, it turns out, that |ma — m(,| = 1 for 
most of the ZL tubes, but there are other cases such as the 
(12, 24) tube. Remember that n = gcd(n, m) and a — \ai\. 

where we have neglected sq in the energy (but not in the 
eigenstate). By inserting k decomposed along the tube 
(T) and around the tube (C): k — k-^ + k^-^ with 

kc = ■]^f^c, one obtains the band structure for nanotube 
labelled by k and ric. Essentially the same tight-binding 
calculation can be done using only the helical and rota- 
tional symmetry as in [jJl and the same result is found, 
when we use k = k • H and Eq. ljA.4|l to convert between 
the quantum numbersiSi (see Fig llOH . 

The Fermi level is at ep = 0, since half of the 
states {2pz orbitals) are filled. By doping and/or a 
gate voltage the Fermi level can be moved about ~ 
ib0.4eVi ^i'''^i'^''i^°i^^i^^ Furthermore note that graphene 
has electron-hole symmetrjiSS, for £p = and therefore 
so does any (n, m) carbon nanotube. 

1. The linearized band structure 

We are only interested in the transport properties of 
nanotubes and therefore expand T(k) around the Fermi 
level = 0, i.e. around the two zeros-'- of T(k), K.; = 
^(^,4), (^ = ±1) and obtain 

T(K, + i^)~ ^(ii^,+^J^^), (A.IO) 

where we have introduced the deviation from by 
Si=k- K,. Note that |T(K, -I- Si)\ ~ ^\Si\ used in 
Eq. (|A.9ll . Furthermore, note that we do not expand T 
around each individual 7^ used, but around Ep = 0, 
since this preserves the electron-hole symmetry of the 
band structure. By inserting K = ^t-^ + ^c-^ into 
the periodic boundary condition Eq. IIA.3p the energy is 
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found to be 



(A.ll) 



2hvo (HtDV ( 



(n + m) + <ri(m — n) 



where D = ^-Vn^+m^+m^ ^j^^ diameter, = is 
the value of the velocity in all metallic tubes. 



2. Unified picture of metallic tubes: armchair- like 
and zigzag-like tubes 



We will now show using the linearized T Eq. ljA.lfl|l 
that all metallic tubes are either zigzag-like or armchair- 
like and define the precise meaning of this. If £ 



1 the (n, m) tube is metallic and has four crossings of 
the Fermi level found from Eq. ijA.llll to be: Two (the 
± in Eq. ifXTlTl l for n'}=+^ and two (±) for 

nl=~^ = (i.e. = 0). This gives the energy and 

eigenstates for the bands crossing the Fermi level: 



— ^hvo^T and 



i{m—n) — <;y/3{n'\-rrL) 
2y/ n'^~{'m-^-\-rnn 

1 



(A.12) 
(A.13) 



where ^ = ±1. This is found by inserting the linearized 
T Eq. l|XT(1l into Eq. and using J?c = 0. By 

doing this, we get sign(.ftT) in a and e oc ±|.^t|, but 
we require continuity of the states (across the T = 
point where Eq. llA.9|l was not valid) and remove the 
sign-function and thereby also the absolute value, i.e. the 
metallic linear bands cross the Fermi level (sp — 0) . Note 
that a and /3 are independent of and thereby k (to 
first order in k), which turns out to be important in the 
Coulomb matrix element. The energy bands cross the 
Fermi level (ep = 0) at .^t = and since k = k ■ = 

.^T + • 1^ the crossing of £p = as a function of k is 
at 



T 

m 



-2Tm f 
gcd(2n-|-m,2m+n)|T| ^ 

27rm f 

gcd(2r!.+m,2m+n)|T| ^ 



(A.14) 



which are either both aX k = (double degenerate, 
nl^^^) or k — i-^^ (non-degenerate) for k in the FBZ, 
— jry < < fry (see 0| for details). Furthermore, 
we have the following connection between the crossing 
of £p = and the crystal angular momentum of the 
bands crossin^^S: If the bands are crossing Cp = at 
k = 0, then the two doubly degenerate crosses have dif- 
ferent nonzero angular momentum tria — ^"^"^ (mod n), 
tTife = 2n^fri ^j^qJ ^j^j -j- jf ^jjg other 



± 



27r 



then both crosses 



hand the crossing is at fc — ---jixi- 
have nXa = mf, = 0. This makes it possible to divide 
all metallic tubes into either armchair-like or zigzag-like 



tubes (see figure I^J with the following bands crossing the 
Fermi level (cp = 0): 



Zigzag-like: 
Armchair-like: 



Chvok, me{ma,mb} (A. 15) 
-mvo{\k\-ko), (m = 0), (A.16) 



where k^ 



^ = ±1, n = ±1 and fee] ^ 



3|T| 



|T|' |T|J 



The translational vector T (and |T|) is different for dif- 
ferent metallic tubes independent of the type. Note that 
the armchair-like bands are in general not connected in 
the way modelled by Eq. ijA.lfip (consider e.g. a (7, 4) 
tube), but since they have the same angular momentum 
m = we connect the bands in this way for convenience. 
For scattering between the bands we will however con- 
sider the bands as four bands as we will see in section 
mil Examples of zigzag-like and armchair-like tubes are 
found in tabled 

For a (real) armchair (n, n) tube the IT index in 
Eq. (|A.16p is the parity in the angular coordinate in cylin- 
drical coordinate a7?i7^ and the states are (j^ — ^ " ^ 

to all orders in k (in the nearest neighbor tight-binding 
approximation). Results similar to the ones obtained 
from the linearized T Eq. (jA.injl can by found by using 
the k • p approximationjiifi however this does not reveal 
the crystal angular momentum. 



APPENDIX B: SCREENING IN THE RPA 
APPROACH INCLUDING THE BAND 
STRUCTURE 

Here we calculate the screened Coulomb potential in 
the random phase approximation (RPA) in order to in- 
clude both static and dynamical screening effects in the 
Coulomb drag, which have been seen to be important 
perviously for bilayer systems i^^'SS'Si 

The Dyson equation for the screened potential in real 
and frequency space is: 

y(ri,r2,c^) =F"(|ri-r2|)+ 

JdrJ dr' y"(|ri - r|)x°(r, r', c^)F(r', rs, cv), (B.l) 

where the non-interacting polarizability is 

X°(r, t, r', t') = -te{t - t'){[p{r, t), p{v' , i')])o, (B.2) 

where /5(r, t) is the density operator in the interaction 
picture and the average (• • • )o is taken for non-interacting 
particles. By writing the density operator by the help of 
a complete set of quantum states {(/3,,(r)} we find the 
polarizability to be 



X°(r,r',c.) 



(B.3) 
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where 0+ is a positive infinitesimal, /^(e) is the Fermi 
function and x°,r;'('^) '^^s introduced. To find the 
Coulomb matrix element we insert the RPA equation 
||RT)| into 

(l'2'|y(ri,r2,c^)|12) = (B.4) 
j dri J dr2ip*i,{ri)(p2,{r2)V{ri,r2,uj)(piiri)ip2{r2) 

and get 

(l'2'|y(ri,r2,c.)|12) = (l'2'|yO(|ri - r2|)|12)+ (B.5) 
Y,xl,'{^){l'v\V'{\ri~r\)\h/){r/2'\V{r',r2,u;)m. 



This equation can be used for any set of quantum 
states and in particular for the metallic states for nan- 
otubes, so rj is the set of indices (i, fc, ^, cr), where 
i = 1,2 is the tube index, a is the spin and remem- 
ber that <; determines the angular momentum m. The 
screened and unscreened matrix elements Eq. Ijlflll and 
Eq. ([Till can now be inserted into Eq. ljB.5p to get the 
screened matrix element. Doing this, we observe that 

5i(fci^i?i,fckiC'i)52(fc2^^2C2:^2^2C2) IS & commou factor, 
which simplifies the result. To simplify further, we use 
that gi and Xi ^ire periodic in the reciprocal lattice 
Gi for subsystem i, g{r],ri') — g*{r]',r]) and introduce 
Qi = k^ — fcj, Anxj = m ■ — tUj and 



>Vi„;2(gi, Ami,g2, Am2,t^) = ^ ^ V{qi + G^,,Ami + n^^Ui^,q2 + Gi^, Am2 + n,^u^^,ri^,ri^), (B.6) 

where ii, 12 are tube indices. Equivalently we introduce Wj^j^ the sum over (without the g's and the 2^ 
factor). So Eq. ljB.5|l becomes 



>V,^,,(gii, Ami,,gi2, Ami2,cj) = 27rL>V°,^((jii, Amii,?^^, Am^J + Xl^^^'^^i + 

Am,,,c^), (B.7) 

which has a matrix structure in the reciprocal lattice and in i and the effective polarization is 



Am, c.) = — ^ ^ x^ik^,, k + q^,', u;)\g,ik,^, ,;k + q, ,')\', 



(B.5 



where ?' is chosen such that m' = m + Am. Note that x° is diagonal in the tube index i, since we do not include 
tunnelling between the tubes. In order to find the screened intershell Coulomb interaction we truncate Eq. IIB.7|l and 
only include the G^^ — and — term in the sum, which gives us a 2 x 2 matrix equation (in i) to find W12, and 
therefore the screened Coulomb matrix element is: 



{k[m[^[ , /c>^^^ |y (ri , r2 , w) IfciTTXi^ , ^32^2^2 



2ttL 



9liki^l^l:k[''[Qg2{k2^2^2^K'^2Q 



Gl,G2 Ui,U2 



V°{k[ - fci + Gi,m[ - ni]^ + n^u^,ri,r2] 
ei2{k[ - ki,m[ ~m^,uj) 



with 



ei2{q, Am, uj) = 1 - XesAl: ^m, uj)V°{q, Am, n, n) 1 - Xeff,2(9, Am, uj)V°{q, Am, r2, r2) 
- Xeff,i(9> Am,u;)x2fF,2(9, Am,u;) x V"{q, Am,ri,r2)V"{q, Am,r2,ri), 



(B.9) 



(B.IO) 



where we have neglected the reciprocal lattice vectors different from zero and therefore used 
^"fei, Am,^,ri,r2)W^2fei> Am,,,g,,,Am,J - V°iq,^ , Am,^ , r2, r2)y vS|2(gn^ Am,, , g^,, Am,J ~ 0. 

If we consider armchair-Hke tubes (only the linear bands from Eq. ljA.16|l l. then all the crystal angular momentum 
is zero and from the g-factor analysis in section UlI Al the interband transition (11 = 1 ^ 11' = — 1 in Eq. IjA.lfiH can 
safely be neglected and for the intraband transition we have 5 ~ 1. Therefore 



°ff,,(g,0,o.) = — ^ J2 X^{kU,k + qU,u;)^xlsAq:U;] 



n=+i 



n=-i 



(B.ll) 



k n=±i 
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and for < q < we find in the long tube limit and for zero temperature (T = 0) 

(fco 



n=+i _ 



(27r)2?i 



2£p 



~ ^ 9 2 2 2 



g - fco + t:^ 



_ ^2 



2^0 



■In 



-0 (q-2k,+ '^) 



9 ^ ofco - -r^ I — In 



— Wgg^ 



and 



(27r)2fi 



\ 2"og ( fco + ^ - 

^" + ? 2"^ T 



fco 



^ — In 



«o U-2fco-& 



(B.12) 



2voq (ifco + ^ - 9 



— VqQ'^ 



2 ^° ^ fiwo . 



2wo9 ( |f| - 



2vo 



In 



^'o 19 



2ep 



w2g2 



(B.13) 



r 



which for small q and w simplifies to the result iniii: 

4voq^ 



(B.14) 



Note that in the static limit the effective polarizability is 
just a constant. The zero temperature approximation^ 
of the polarizability is good as long as T is much smaller 
than Tp, which is often the case for nanotubes (Tp ~ 
lOOOK). Including finite temperature in the polarizability 
could give a plasmon enhanced drag as previously found 
for bilayer systems^SiS at T ~ 0.5Tp. 

For zigzag-like tubes the effective polarizability can 
be found in the same way, but for the linear bands 
crossing the Fermi level (sp = 0) we can - in con- 
trast to the armchair-like case - have both Am = and 
Am = ±(ma - mb). 



The unscreened Coulomb interaction V^{q, Am, ri,rj) 
can be found from the Poisson equation by Fourier trans- 
forming in the cylindrical coordinate and in the coordi- 
nate along the tube, i.e 



80 



V°{q,Am,n,rj) = —lAm{qr^)KAmiqrJ) n < r,, (B.15) 
Co 

where lAm(a;) (KAm(a;)) is the modified Bessel's func- 
tions of the first (second) kind of order Am and eo is the 
vacuum permittivity. Note that the small q limit is (log- 
arithmic) divergent only for the potential with Am = 0. 

So we have all the ingredients in the screened Coulomb 
matrix element between different shells using the tight- 
binding states of the carbon nanotubes, which is used to 
model the Coulomb drag between the shells. 

*E-mail: lunan@fys.ku.dk 
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